Q^ '■ Disorder Effects in CA-Models for Traffic Flow 

Q ' W. Knospe^, L. Santen^, A. Schadschneider^, M. Schreckenberg-'^ 

(D ' 

Q ■ ^ Theoretische Physik FBIO 

Gerhard-Mercator Universitat Duisburg 
D-47048 Duisburg, Germany 



^ 



^ Institut fiir Theoretische Physik 
f^ ■ Universitat zu Koln 

O ■ D-50937 Koln, Germany 



Cw . Abstract. We investigate the effect of quenched disorder in the Nagel- 

f/3 ' Schreckenberg model of traffic ffow. Spatial inhomogenities, i.e. lattice 

sites where the braking probability is enlarged, are considered as well as 
Cu ■ particle disorder, i.e. cars of a different maximum velocity. Both types of 

S ' disorder lead to segregated states. 
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1 Introduction 



The formation of traffic jams is one of the fundamental problems of traffic flow 
theory. Traffic jams can form spontanously as well as due to hindrances, e.g. 
K^ ■ road works or slow cars. Although these hindrances often cover only a small 

J^ , part of the road, they can cause large jams, or, in a more physical language, one 

observes macroscopic effects due to local defects in the system. This behaviour 
p^ \ is one of the characteristic properties of nonequilibrium systems. In the context 

of driven lattice gases one can distinguish in principle two types of defects: 1) 
C""^ ! lattice defects, e. g. sites where the mobility of particles is reduced, and 2) particle 

^*^ ■ defects, for example 'slow' particles. 

For the simplest model for a driven lattice gas in one dimension, the asym- 
metric exclusions process (ASEP)[n|, it has been shown that both types of defects 
can cause phase transitions. Slow particles, i.e. particles with a reduced hopping 
'^ ■ rate, determine completely the flow at small densities. Moreover, if the distri- 

M , bution of the hopping rates fulfills certain conditions, one can show analytically 

O ' that the variance of the distance distribution has a logarithmic divergence at the 

transition density (|,|,|,D • 

Also lattice defects in the ASEP have been investigated extensively [|6|J7|,|^,p|, [lO|jll|| . 
Implementing lattice defects one can observe three different phases: A high and 
^ ' a low density phase, where the average flow of the homogeneous system is re- 

5^ , covered, and a segregated phase at intermediate densities, where the flow takes 

a constant value. Despite extensive effort, exact analytical results exist only for 
a special type of update |13] , but approximative descriptions are in reasonable 
agreement with simulation results. 



In this work we generalize the resuhs of the ASEP to the Nagel-Schreckenberg 
(NaSch) model of traffic flow (for a detailed explanation of the model see [13,1^). 
Compared with the ASEP, the particles (cars) can hop more than one site in 
a single update step and the update rules are applied in parallel to all cars. 
Therefore slow cars can be considered in two ways. First, one can think of cars 
with a smaller maximum velocity and second, of cars with an enlarged braking 
probability. This has been done very recently and analogous results to the v^ax = 
1 case have been found fl^ . In contrast to fl^ we discuss the case of two different 
maximum velocities in the third chapter. In the next chapter we show simulation 
results for the NaSch model with defect sites, implemented as sites where the 
braking probability of cars is higher compared to the rest of the lattice. An 
alternative choice has been used in |lq| , where a 'speed limit' in a part of the 
lattice has been considered. 



2 Defect Sites 

As already mentioned above, in this chapter we show simulation results of the 
NaSch model on a lattice with defect sites. Fig. Ill shows the fundamental diagram 
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Fig. 1. Fundamental diagram of the NaSch model with a single defect site. The 
model parameters are given by Vmax = 2 (maximum velocity), p = 0.25 (braking 
probabihty), pd = 0.75 (braking probability at the defect site) and L — 3200 
(system size). 



of a system with a single defect site. Obviously we can distinguish three different 
phases depending on the density. In the high and low density phases the average 
flow of the defect systems takes the same value as in the homogeneous system. 
For intermediate densities the flow is constant and limited by the capacity of the 
defect site. 

This behaviour of the average flow can be explained looking at the density 
proflle. In the high and low density phase only local deviations from a constant 
profile can be observed, but at intermediate densities one observes a separa- 
tion into macroscopic high and low density regions. Changing the global density 
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Fig. 2. Density profiles in the three phases. In the high (red, p = 0.60) and low 
(blue, p = 0.10) density phases only local inhomogenities occur near the defect 
site (the defect site is located at Xd = 1600 = L/2), but for intermediate densities 
one observes phase separation (green, p = 0.30). We used the same parameters 
as in Fig. |l|. 



within the segregated phase, the bulk densities in the high (ph) and low density 
region (pi) remain constant, only the length of the high and low density region 
changes. Consequently the average flow is constant in the segregated phase, be- 
cause the average density in the vicinity of the defect site does not depend on 
the global density. 

Near the average position of the shock the density profile decays exponentially 
from ph to pi. Therefore one can introduce a length scale ^, which corresponds 
to the magnitude of the fluctuations of the shock position. For Vmax > 1 we 
found numerically that ^ ~ VL for all densities in the segregated phase we took 
into account. This scaling behaviour is already known for the ASEP at p 7^ 0.5. 
The modified scaling behaviour ^ ^ L^^^ at p = 0.5 was not found for Vmax > 1- 
This confirms the picture that the reduction of fluctuations is a consequence of 
the particle- hole symmetry ||^. 

A good estimate for the plateau value of the flow for the case of Vmax = 1 
can be obtained using the assumption that the system is separated into two 
regions of constant density, where the results from the homogeneous system can 



be used. Using an argumentation similar to 1 10 we obtain the plateau value of 
the average flow: 
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The bulk value of the density in the high (low) density region the density is 
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treatment of the defect system is in a resonable agreement with the simulation 



results for large system sizes (for a more detailed discussion see llq] ) . In principle 
one could treat higher velocities in the same way, but unfortunately no exact 
analytical description of the homogeneous system has been found by so far. 



3 Particle defects 

In order to study the effect of slow cars on the throughput of single lane traffic, 
we discuss a system which consists of fast cars with maximum velocity Vf = 3 
and one slow car with v^ = 2. 




Fig. 3. Fundamental diagram of a system with one slow car (dots) compared 
with fundamental diagrams of the corresponding homogeneous systems. The 
maximum velocity of the slow car is given by z;s = 2 and of the other cars by 
Vf — 3. The braking probability of all cars is p = 0.5. 



In Fig. ^ the fundamental diagram of this system is compared with the fun- 
damental diagram of the analogous homogeneous systems. Obviously for low 
densities the flow is given by Ji{p) — p{vs — p) in a agreement with the ho- 
mogeneous system with Vmax = 2. Compared to the homogeneous system one 
observes the linear density dependence also for larger values of the average den- 
sity. The average flow depends linearly on the global density if Ji{p) < Jv;{p) 
holds, where Jv; (p) denotes the stationary flow of the homogeneous system with 
fast cars. Deviations from this form that are observable near the intersection 
point of Ji[p) and Jv^{p), are due to the finite size effects. 

The reason for the stability of the free flow regime in the mixed system is the 
gap in front of the slow car in the stationary state. In Fig. |4| the density profile 
relative to the position of the slow car is shown. Obviously a large jam behind 
the slow car appears and the density in front of the slow car vanishes. Therefore 
the slow car can move with its free flow velocity Vg — p. This velocity determines 
the velocity of the whole jam and therefore the average velocity of the cars is 
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Fig. 4. Density profile relative to the position of the slow car. The data are 
obtained for the following set of parameters: Vs = 1, Vf = 2, p = 0.5 and, 
p = 0.1. 



simply given by Ji{p). This is one possible simple scenario for the occurance of 
'moving' jams that are well known from measurements. 

In addition to the macroscopic form, two details of the density profile have 
to be discussed. Behind the slow car one observes an oscillating amplitude of the 
density profile. These oscillations are well know from the spatial correlation func- 
tion of the homogeneous system for small p and densities near pc — l/{vmax + 1) 
[ pjj| . Obviously the slow car synchronizes the motion of the faster cars. Further- 
more we want to discuss the decay of the density profile at the end of the jam. 
For finite values of the braking probability the jamlength is fluctuating. These 
fluctuations are responsible for the finite size effects discussed above, because 
temporarily the gap in front of the slow car vanishes due to fluctuations of the 
jam length. Measurements of the jamlength show that the fluctuations scale ac- 
cording to ALj ^ \fTj. Due to the subextensive scaling of the fluctations the 
segregated states are observable up to p = pj if L — > oo, where pj denotes the 
bulk density of the jam. 



4 Summary 

The investigation of local defects in the NaSch model has shown that these de- 
fects can change the macroscopic properties of the model. In the case of lattice 
defects three phases occur analogous to the ASEP (i.e. Vraax — 1) with continu- 
ous time update. The segregated phase also exists for higher maximum velocities, 
but some details of the density profile are changed. In the low density regime 
one observes an oscillating profile as a consequence of the parallel update. The 
anomalous scaling of the fluctuations of the shock position found for the ASEP 
at yO = 0.5 is absent for Vmax > 1, for all parameter combinations that have 
been taken into account. This result confirms the picture that the particle-hole 
symmetry is essential for the reduction of fluctuations. 



Particle defects also produce phase separated stationary states at low densi- 
ties. These states consist of a large jam behind the slowest vehicle and a large 
gap in front of the slowest car. The velocity of the 'moving' jam is completely 
determined by the free flow velocity of the slow car. The formation of large clus- 
ters due to slow cars can also be observed in two-lane traffic, but due to lane 
changing the lifetime of the large clusters is finite (at least for low densities). 
Nevertheless it has been shown, that already a small concentration of slow cars 
changes leads to a drastic reduction of the average flow at low densities [pOl . 

Our results show that the system poperties can change completely, if local 
disorder is taken into account. These results might also be important for an 
understanding of the behaviour of real traffic which are often determined by 
'imperfections'. Such imperfections are not only the defects considered here, but 
also other deviations from an ideal system, e.g. ramps or a finite system size. 
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